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ABSTRACT 

We discuss the rest-frame UV emission between 5 < z < 7 from the MareNostrum 
High-z Universe, a SPH simulation done with more than 2 bilhon particles. Cosmologi- 
cal simulations of galaxy formation generally overpredict the UV restframe luminosity 
function at high redshift, both at the bright and faint ends. In this Letter we explore 
a dust attenuation model where a larger extinction is applied to star populations 
younger than a given age, mimicking the effect of a clumpy interstellar medium. We 
show that this scenario fits reasonably well both the UV luminosity functions and the 
UV-continuum slopes derived from observations. The model assumes a large obscura- 
tion for stars younger than 25 Myr from the gas clouds where they should be embedded 
at their formation time. We find that the optical depth in these clouds should be be- 
tween 30 and 100 times larger than the mean optical depth for the homogeneous part 
of the interstellar medium. These values are one order of magnitude larger than those 
estimated in local galaxies. Therefore, we conclude that ACDM predictions for the 
high-z UV emission can accommodate the current observations if we consider a dust 
extinction model based on the assumption of a clumpy environment at high redshift. 



One of the most relevant issues in cosmology is under- 
standing the evolution of galaxy populations. Thanks to 
recently developed techniques and improved instrumenta- 
tion, observations of high redshift gala xies (z > 5) be- 
gan to be possible during the l ast years ( Iwata et al.l 20071: 



Bouwcns ct al. 2007, 2008; Mc Lure et all 12003 : IStark et al 
2009; O uchi et al.i 2009: Bou wens et a l' 2009a b- IOesch et al 



2009. : .Bunker et al..,2009. ; ,McLure et al...2009 ') . 

The fundamental quantity used in the exploration of 
these galaxy populations is the galaxy luminosity function 
(LF). In particular, the LF in the rest-frame UV provides 
valuable information on key physical processes. The UV 
output of a galaxy imposes constraints on the instanta- 
neous formation rate and amount of young stars in a galaxy. 
This kind of census provides information on the role of 

lalactic star populat ions in the reionization of the Universe 

Stanwav et al.ll2003l '). 

Observationally, the method of choice to gather a large 
sample of galaxies at high redshift is the drop-out technique. 
This technique detects the spectral discontinuity at 912 A 
(Lyman break) due to the absorption of UV radiation by in- 



tergalacti c hydrogen gas usin g multi- wavelength broad band 
imaging (I Steidel et al.l 1 19961 ) . The galaxies discovered us- 
ing this method receive the name of Lyman Break Galaxies 
(LBGs). 

Motivated by the observational results, recent theoret- 
ical studies of LB Gs beyond z > 5 have been developed 
iNight et al.l 1200^ ). These theoretical studies, usually re- 
quire computational techniques to produce a mock galaxy 
population which can be directly compared with the ob- 
served one. In the case of the rest-frame UV LF, two impor- 
tant elements are required in the model: the stellar popu- 
lation (masses, ages and metallicities) and the correspond- 
ing gas content in order to derive the extinction proper- 
ties. Modelling the intrinsic UV emission through stellar 
population synthesis models is now a standard procedure 
(Maraston 1998). On the other hand, having a proper esti- 
mation of the extinction is a more uncertain process. 

The estimation of th e extinction c an assume the form 
of a simple Calzetti law (iCalzetti et al.l i2000'l , where usually 
the UV magnitudes of each galaxy are reduced by the same 
amount reg ardless of the ph ysical properties of the gas in 
the galaxy l|Night et al.ll200^ ). If the galaxy model provides 
information on the gas and its metallicity contents, it would 
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be desirable to implement a more physical model for the dust 
obscuration and try to put some constraints on the amount 
and distribution (i.e. dumpiness) of the dust on galactic 
scales (|lnouc 2005). 

In this Letter we present the first results from the 
MareNostrum High-z Universe. The simulation follows self- 
consistently the dynamical evolution of more than 2 billion 
dark matter and gas particles and includes a prescription for 
star formation, supernovae feedback and UV background. 
Here, we will present the numerical estimates for the con- 
tinuum rest-frame UV emission from the stellar component, 
and a dust attenuation model. Our main objective is to com- 
pare the numerical predictions with the observed evolution 
of the UV LF between redshifts 2 = 5 and z = 7. 



analysis in AHF. The thermal energy of gas particles is taken 
into account during the calculation of the binding energy. 

All objects with more than 1000 particles, dark matter, 
gas and stars combined, are used in our present analyses. 
We assume a galaxy is resolved if the object contains 200 
or more star particles, which corresponds to objects with > 
400 particles of gas. This ensures a proper estimation of the 
average gas column densities in our numerical galaxies. We 
then build the full Spectral Energy Distribution (SED) and 
compute the dust attenuation using the information from 
the stars and gas particles of each galaxy. In the following 
section we explain the details of these procedures. 



2 SIMULATION AND GALAXY FINDING 

The MareNostrum High-z Universe simulatiorQ follows the 
non linear evolution of structures in baryons (gas and stars) 
and dark matter, starting from z = 60 within a cube of 
50/i"^Mpc comoving on a side. The c osmological parame - 
ters used correspond to WMAPl data (jSpergel et aklbOOSl ) 
and are = 0.3, fib = 0.045, = 0.7, as = 0.9, a Hubble 
parameter h = 0.7, and a spectral index n = 1. The initial 
density field has been sampled by 1024^ dark matter parti- 
cles with a mass of moM = 8.2 x IO'^/i'^Mq and 1024^ SPH 
gas particles with a mass of mgas = 1.4 X lO^/i^^Mg. The 
simulation has b een performed using the TREEPM-I-SPH 
code GADGET-2 jSpringell l2005h . Radiative and Compton 
cooling processes for an optically thin primordial plasma 
of Helium and Hydrogen are included. We assumed photo- 
ionisatio n by an external spatially uniform UV-background 
adopting iHaardt fc Madciul (|l996l ). The physics of star for- 
mation is treated by means of a sub-resolution model in 
which the gas of the interstellar medium (ISM) is descr ibed 
as a multiphase medium of hot a nd cold gas ( Yepes et al.l 
Il997l : ISpringel fc Hernguistl l2003h . Stars can be formed 
in regions that are sufficiently dense and cold. We con- 
sider the thermal feedback of supernovae as well as the 
effects of stellar win ds following the model described in 
ISpringel fc Hernguistl (|2003) The gravitational smoothing 
scale was set to 2 /i~^kpc in comoving coordinates. The sim- 
ulation was run in the MareNostrum supercomputer using 
800 processors simultaneously. This simulation is intended 
to study the early phases of galaxy formation. It is currently 
at z=4.55, after spending more than 4.5 million cpu hours. 

We identify the objects in the simulations using the 
AMIGA Halo Finder (AHpj^ which identifies both halos and 
subhalos. AHF is an improvement of the HHP halo finder 
iGill et al.l |2004| ). which locates local overdensities in an 
adaptively smoothed density field as prospective halo cen- 
ters. AHF is described in detail in iKnoUmann fc Knebd 
l|2009ll . 

In AHF the local potential minima are computed for each 
of these density peaks and the gravitationally bound par- 
ticles are determined. Only peaks with at least 20 bound 
particles are considered as haloes and retained for further 



^ |http : //astro . f t . uam . es/marenostrum' 

It is a MPI-|-OponMP hybrid halo finder to be downloaded 
freely from htt p: //www.popia.f t .nam. es/AMIGA| 



3 SPECTRAL MODELLING 

The photometric properties of the galaxies are calculated 
employing the st e llar p opulation synthesis model STARDUST 
JPevriendt et al' '19991). using the methods described in 
iHatton et al. (2003). We adopt a Salpeter Initial Mass Func- 
tion (IMF). The SEDs are built using the AHF catalogs. Nev- 
ertheless, we have checked that all the results quoted in this 
Letter remain valid if we identify galaxies by the Friend-of- 
Friends algorithm that run over the star particles only. 

The dust attenuation model parametrizes both the ex- 
tinction in a homogeneous interstellar medium (ISM) and 
the molecular clouds around young s tars, following the phys- 
ical model of ICharlot fc Fail l|2000l '). The attenuation from 
dust in the homogeneous ISM assumes later a slab geometry, 
while the additional attenuation for young stars assumes a 
spherical symmetry. 

We describe first the optical depth for the homogeneous 
interstellar medium, denoted by ■ We take the mean per- 
pendicular optical depth of a galactic disc at wavelength A 
to be 



H 



{Nh) 



2.1 X lO^^atoms cm 



(1) 



where A\/Av is the extinction curve from iMathis et al.l 
(|l983h . Zg is the gas metallicity, {Nh) is the mean atomic 
hydrogen column density and r7=(l + 2)~"isa factor that 
takes into account the evolution of the dust to gas ratio 
at different redshifts, with a > from the a vailable con - 
straints based on simplified theo retical models lllnouell2003l ) 
and observations around z ~ 3 (|Reddv et al.ll2006l ). Given 
the uncertainties in this value, we will treat a as a free pa- 
rameter. 

The mean H column density is calculated as 



{Nh) = Xt_ 



Ma 



-atoms cm 



(2) 



where Xh ~ 0.75 is the universal mass fraction of Hy- 
drogen, Afg is the mass in gas, rg is the radius of the galaxy 
and nip is the proton mass. The extinction curve depends 
on the gas metallicity Zg and is based on an interpolation 
between the solar neighborhood and the Large and Small 
Magellanic Clouds (r = 1.35 for A < 2000A and r = 1.6 for 
A > 2000A). 

The radius, stellar and gas masses for each galaxy are 
taken from the AHF catalogs. Computing {Nh) from the 
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galaxy catalogs turns out to give similar results, on aver- 
age, than integrating the 3D gas distribution of the galaxies 
using the appropriate SPH kernel. The difference between 
the average Hydrogen column densities calculated in these 
two different ways axe less than 15% regardless of the size 
of the galaxy, provided that the gas distribution is sampled 
with more than 200 particles. 

In addition to the foreground, homogeneous ISM extinc- 
tion, we also model, in a simple manner, the attenuation of 
young stars that are embedded in their nest clouds. Stars 
younger than a given age, tc, are subject to an additional 
attenuation with mean perpendicular optical depth 

(3) 

where p is the fraction of the total optical depth seen by 
the young s tars that is found in the homogeneous ISM. 
iKong et al.l (|2004l ) in a study of 115 nearby galaxies, span- 
ning a wide range of star formation activities, find that ^ 
is typically ~ 1/3 with a wide range of scatter between 0.1 
and 0.6. Given the possible uncertainties in extrapolating 
the typical values of to high redshift, we will treat this 
quantity as a free parameter. 



4 ANALYSIS AND RESULTS 

In the calculation of the luminosity functions we have made 
first an additional correction motivated from the different 
abundance of dark matter halos between t he cosmology 
used in the simulation (WMAPl, erg = 0.90 (|Spergel et all 
and the values of the cosmological parameters esti- 
mated from t he latest CMB results (WMAP5, erg = 0.796 
iDunklev et al . 2009)). At a given redshift and mass range, 
the WMAPl universe has a larger halo number density 
than in the WMAP5 one. We estimate the influence of 
a different cosmology u sing the Sheth-Tormen formalism 
l|Sheth fc TormenI 1 19991) . The ratio between the WMAP5 
and WMAPl halo mass functions is used as a weight factor 
for each galaxy (depending on the mass of the host halo) 
in the construction of the LF. This procedure is justified 
because of the strong correlation shown between the UV 
luminosity and the halo mass. 

We have tested this correction by running two addi- 
tional simulations with exactly the same mass and force res- 
olutions but a corresponding smaller volume to speed up the 
calculations (i.e 2 x 256"^ particles in a 12.5 h-^Mpc comov- 
ing boxsize). They were performed with same initial condi- 
tions using the WMAPl and WMAP5 cosmological param- 
eters. Though not shown here, we do reproduce the results 
of the WMAP5 run when using the WMAPl run with the 
aforementioned correction. 

The results after correcting for the halo abundance are 
shown in Figure [1] The black squares represent the results 
without any dust extinction. All the rest-frame UV LFs, 
without dust extinction, strongly over-predict (up to 2 or- 
ders of magnitude) the observational estimates in all the 
redshift range, in fa irly agreement wi th results from other 
ACDM simulations (|Night et al.ll2006l ). 

This discrepancy cannot be attributed to the selection 
of the Salpeter IMF to compute the SEDs. Other IMFs, 



such as Chabrier (|Chabrieij|2603l ) or Kroupa l|Kroupall200l ) 
have lower mass to luminosity ratios in the UV (at most 
~ 1.5) with resp ect to a Salpeter IMF, for a youn g starburst 
of a given age. (|Chabried |2003 : iBruzual fc Chariot 200i). 
These IMFs would then increase the discrepancy between 
our simulation and observations. 

One could argue that the excess of galaxies in our sim- 
ulation could be related with some of the astrophysical pro- 
cesses that shape t he luminosity function at low redshift 
(| Benson et al.l [20031 ) . In particular, an improved modelling 
of supernovae and quasar feedbacks can modify the lumi- 
nosity function. This scenario cannot be ruled out on the 
basis of the present simulation, and certainly deserves fur- 
ther study. Nevertheless, in our case, we want to show that 
a simple, physically motivated extinction model can explain 
the different behaviour of the simulated and observed UV 
LF. 

If we apply only the extinction due to the homogeneous 
ISM, parametrized by = 1.0 (for any value of a) then only 
the most massive galaxies are effectively extinguished. If we 
modify the values of a in such a way as to match the bright 
end of the observed LF, we still find an excess of simulated 
galaxies at the faint end. To modify the luminosities at the 
faint end, specially at z = 5 and z = 6, one needs to apply a 
larger extinction on the young stars. The correction on the 
population of young stars turns out to be independent on 
the total mass of the galaxy, modifying both the bright and 
faint ends of the LFs. 

We explore the parameter space in order to find the best 
set of parameters (a, /x, tc) in our extinction model which 
provides a good fit to the observational data. In Figure [T] 
the red hexagons show the results for the best fit parame- 
ters: a = 1.5, fi — 0.01 — 0.03 and tc — 25Myr. It means 
that stars younger than tc — 25Myr are almost completely 
extinguished. Thus, the optical depth in the clouds embed- 
ding the young stars must be between 30 (/i = 0.03) and 100 
(/i — 0.01) times larger than the optical depth in the ISM. 
For compar i son, th e largest values found in local galaxies by 
iKong et al.l l|2004l ) are ~ 10 times the optical depth in the 
ISM. Furthermore, we do not find any strong degeneracy 
between the parameters a and fi. 

We consider now the colors of the galaxies. One of the 
simplest ways to observationally estimate the dust attenua- 
tion is through the UV slope (3 {fx oc A''). In Figure [5] we 
show the values for f3 in the simulate d galaxies and compa red 
to recent observational estimations l|Bouwens et al.ir2009l lbl). 
We find that that the extinguished colors are consistent with 
the constraints worked out from the observations. There is a 
minor discrepancy at the faintest magnitudes (Mas ~ -19) 
which is not highly significant if we consider that the error 
bars in the Figure [2] represent a 1-a scatter. Remarkably, 
some extinguished simulated galaxies at z ~ 7 (not shown 
in the Figure) still show very blue slopes /3 ~ 3 as seen in 
the observations. 



5 CONCLUSIONS 

We have used the MareNostrum High-z Universe, currently 
the largest SPH simulation of galaxy formation, with more 
than 2 billion particles, to study the evolution of the rest- 
frame UV luminosity function between 5 < z < 7. We find 
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z=5.0/WMAP5 corrected/tc = 25 Myr 



z=6.0/WMAP5 corrected/tc = 25 Myr 



z=7.0/WMAP5 corrected/tc = 25 Myr 
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Figure 1. Luminosity functions at redshifts z 5, 6 and 7. The error bars are poissonian. The horizontal hne shows the limits of one 
object per unit magnitude in the simulation volume, for the used bin size AMab = 0.50. The vertical line indicates the approximate limit 
for well resolved objects following a criteria of more than 200 star particles. This limit corresponds to similar UV magnitudes at every 
redshift. The simulation data has been already corrected for the different abundance of dark matter halos between the WMAPl and 
WMAP5 cosmologies. The open symbols represent the observations. The filled symbols represent different extinction prescriptions on 
the simulated galaxies: no extinction (black squares), homogeneous ISM plus an intermediate extinction on young stars (blue triangles) 
and a extreme extinction on young stars (red hexagons). In the model with extreme extinction on young stars, all the stars younger than 
25 Myr are almost completely extinguished. 
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Mab(UV) Mab(UV) Mab(UV) 

Figure 2. The slope, /3, of the UV-continuum versus absolute UV magnitudes of simulated galaxies at 2 ~ 5, 6 and 7. The blue points 
correspond to the /3 determinations without extinction, red points correspond to the estimation with the extinction model consistent 
with the UV luminosity functions. The open symbols represent the estimated values of f3 from HST observations. In both cases the 
vertical error bars represent a l-cr scatter. The brightest galaxies in the simulation follows the same color trend as in observations. Some 
discrepancy is still present at the faint end though. 



that the UV LF from the simulation over-predicts system- 
atically the observations. Using a Chabrier, Kroupa or Top 
Heavy IMF, instead of the Salpeter IMF we have used here, 
will increase the disagreement. 

We further correct the luminosity of each galaxy by dust 
attenuation. If we correct only by the attenuation produced 
by an homogeneous ISM, with an optical depth calculated 
using Equation [1] we find that the results are unsatisfactory. 
In that kind of extinction model, the dust optical depth 
turns out to be dependent on the mass of the galaxy, making 



the bright end of the LF more extinguished than the faint 
end, the latter remaining basically unchanged. 

We use then a simp l e but physically motivated model 
based on lCharlot Sz Fail (HqoO), where the extinction is de- 
scribed by a homogeneous ISM component plus a clumpy 
component affecting the younger stars. We find that the 
additional extinction on the younger stars is necessary to 
reproduce the overall normalization and shape of the lumi- 
nosity functions between 5 < z < 7. The final modification 
on the UV luminosities is independent on the total mass of 



High-z UV: hints for a clumpy ISM? 5 



the galaxy, meaning that the faint end of the LF can be 
modified as weU. 

We find that all the stars younger than 25 Myr must be 
obscured almost completely in order to get a good agreement 
to the observed evolution of the LFs. This means that the 
optical depth in the clouds embedding the young stars must 
be between 30 and 100 times larger than the optical depth 
from the homogeneous ISM. Compared to observations of lo- 
cal galaxies this proportion between the two optical depths 
is one order of magnitude larger than the highest estimates. 
The agreement of our model to available observational con- 
straints hints for a clumpy ISM at high redshift. 

The extinction model we propose is consistent with 
the observational constraints of the UV slope between 5 < 
z < 7. We even find some galaxi es with very blue UV 
slopes as seen in galaxies at 2 ~ 7 l|Bouwens et al ] |2009bt ). 
The agreement is less satisfactory at lower UV magnitudes 
Mab ~ —19. The extinction model we apply to the very 
young stars is thus not visible as a strong change in the 
UV slopes, but could be observationally detected in the 
rest-frame infrared with instruments like HERSCHEL and 
ALMA. A complementary view of the results presented so 
far will be provided in an upcoming publication on the phys- 
ical nature of Lyman-a Emitters in the MareNostrum High-z 
Universe. 
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